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The magnetic properties and the effects of interparticle interaction on it have been studied in 
nanoparticles of half doped Pro.sSro.sMnOs. Three samples consisting of nanoparticles of different 
average particle sizes are synthesized to render the variation in interparticle interaction. Though all 
the samples crystallize in the same structure to that of their bulk compound, the low temperature 
ferromagnetic-antiferromagnetic transition, which is present in bulk compound, is not evident in the 
nanoparticles. Linear as well as nonlinear ac susceptibility coupled with dc magnetic measurements 
have shown the superparamagnetic behavior of these nanoparticles where the blocking temperature 
increases with the increasing particle size. Presence of interparticle interaction is confirmed from 
the temperature variation of coercive field and the analysis of frequency dependent ac susceptibility. 
We have identified the nature of this interaction to be of dipolar type, and show that its strength de- 
creases with the increasing particle size. The effect of this dipolar interaction on magnetic properties 
is intriguing as the compounds exhibit crossover from de Almeida- Thouless to Gabay- Toulouse like 
critical lines on field-temperature plane above their respective interaction field. In agreement with 
theoretical prediction, we infer that this crossover is induced by the unidirectional anisotropy arising 
from interparticle interaction, and this is confirmed from the presence of exchange bias phenomenon. 

PACS numbers: 75.75.-c, 75.47.Lx, 75.40.Gb, 75.20.-g 



I. INTRODUCTION 

The complex coupling between spin, charge, orbital 
and lattice in perovskite manganites around half doping 
gives rise to various exotic electro-magnetic properties^— 
The of critical balance between the competing antifcrro- 
magnetic (AF) - insulating (I) and ferromagnetic (FM) - 
metallic (M) orders is found to be rather susceptible to 
small perturbations which can causing drastic effects^^ 
In the recent time, it is shown that reduction in the par- 
ticle size or dimensionality have rather significant effects 
on the physical properties. The charge ordering (CO), or- 
bital ordering (00) and low temperature (T) AF states 
are largely modified in the thin films of Ndo.sSrrj.sMnOs 
and Pro.5Sro.5Mn03<^ii In the parallel developments, it is 
demonstrated that in the nanoparticles, the notable bulk 
properties like: CO, 00, low temperature AF states are 
largely suppressed or weakened, and the system emerges 
to FM stated— However, it was suggested by Neel long 
back that small or finite-size particles of AF material 
should exhibit supcrparamagnctism (SPM) or weak FM 
properties J£ 

The physical properties in nanoparticles are mainly 
modified due to (i) the finite size effect, and (ii) the 
surface disorder effect which is caused by defects, bro- 
ken exchange bonds, fluctuations in number /separation 
of neighboring atoms at the surface. Usually, FM par- 
ticles with reduced size behave as a single domain en- 
tity exhibiting SPM behavior where the thermal fluctu- 
ation overcome the anisotropy energy of individual par- 
ticle causing the magnetization to flip from one to an- 
other easy direction. The situation becomes complicated 
if there exists interparticle interactions which induces col- 
lective behavior, and may modify the magnetic proper- 



ties drastically. For instance, it is observed for ferrofluid 
that low-T dynamics changes from SPM to spin glass 
(SG) like behavior with just increasing interaction among 
the particles*^ Therefore, it is necessary to identify the 
nature as well as the extent of this interaction to char- 
acterize the magnetic properties of the collection or con- 
glomerate of nanoparticles. Commonly, the inter-particle 
interaction is achieved or modulated with the variation 
of volume concentration of particles, however, the size 
variation of particles offer another avenue in this regard 
which does not require the intervening materials unlike 
the former case, and is uncommon in literature. 

Here we present a detailed study of the magnetic 
properties and their evolution with the inter-particle 
interaction in nanoparticles of half doped manganitc, 
Pr .5Sr .5MnO 3 (PSMO). In this study, we have modu- 
lated the inter-particle interaction by means of modifying 
the size of the particles. Nanoparticles of three different 
average sizes (in the range of 15-27 nm) are synthesized 
which have identical room temperature crystal symme- 
try as their bulk counterpart. However, contrary to the 
bulk PSMO, these nanoparticles do not show the low 
temperature FM to AF transitionj 14 i 15 Using linear and 
non-linear ac susceptibility coupled with dc magnetiza- 
tion, we have confirm that these nanoparticles show SPM 
behavior where the blocking temperature (Tb) increases 
with the particle size. Analysis of frequency dependent 
peak in ac susceptibility and temperature variation of 
coercive field indicate the presence of interparticle inter- 
action. The nature of the interaction is found to be of 
dipolar type, and its strength decreases with increasing 
particle size. It is rather significant that with increasing 
the field (H), all the nanoparticles exhibit a crossover in 
critical lines from the de Almeida-Thouless to the Gabay- 
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Toulouse like behavior around the respective dipolar in- 
teraction field. Following the theoretical prediction, we 
believe that this crossover is induced by unidirectional 
anisotropy which, in the present systems arises from the 
intcrparticle interaction and is confirmed from the pres- 
ence of exchange bias effect. 



II. EXPERIMENTAL DETAILS 

Nanocrystalline samples of PSMO are prepared, us- 
ing Pr 6 On, SrC0 3 and Mn(CH 3 COO) 2 with purity of 
99.99% or better, by chemical (pyrophoric) method^ 
First, the proportional amount of ingredients are dis- 
solved in dilute HNO3. Tricthanolamine (TEA) is mixed 
with this aqueous solution of the ingredients in the metal 
ions to TEA ratio as (Pr,Sr):Mn:TEA =1:1:4, and evapo- 
rated around 200°C. At the end, it yields black precursor 
powders which are heated at 600, 650 and 700° C for 5 
hours to get the powders of different nanosize samples 
which are designated as N600, N650 and N700, respec- 
tively. Room temperature x-ray diffraction (XRD) pat- 
terns are collected with Rigaku Dmax 300 diffractomctcr 
attached to 18 kW rotating anode (Cu) source. Rietveld 
refincmenfeil of XRD patterns show that samples are in 
single phase, and they crystallize in tetragonal structure 
with 14/ mem symmetry similar to their bulk counterpart. 
Fig. la presents XRD data along with the Rietveld fit- 
ting for the N650 sample. Average particle size (Dxrd) 
has been estimated from the XRD plot using Scherrcr 
formula*^ The calculated particle sizes are given in Ta- 
ble 1. Structural parameters found from the Rietveld 
refinements are also given in Table 1. There is no major 
structural modification in nanometric samples compared 
to the bulkJ^ Extreme right column of Table 1 presents 
the percentage change in the related parameters for all 
the samples, showing that no major structural modifica- 
tion occurs with the size variation of particles. 

Samples are also characterized through transmis- 
sion electron microscope (TEM) (model: TECNAI G2- 
20FEI). Fig. lb shows TEM bright field image of N650 
sample. Histogram for particle size distribution is ob- 
tained after analyzing several TEM images. Such his- 
togram for N650 sample, presented in Fig. lc, shows that 
particles are polydisperse in nature. Mean particle size 
{Dtem) for all the samples have been obtained from the 
Gaussian fitting of figure similar to Fig. lc and values are 
given in Table 1. The values of Dxrd and Dtem match 
reasonably well except for N700 sample. This mismatch 
in sizes for N700 may be arising from conglomeration of 
particles because of higher temperature processing. Se- 
lected area electron diffraction (SAED) graph has been 
given in Fig. Id for N650 sample. The figure shows 
well distinct spotty concentric ring indicating crystalline 
nanoparticles. The similar figure for N600 sample is given 
in Ref. [l9|. The oxygen stoichiometry or Mn 3+ /Mn 4+ 
ratio of 1:1 is ensured from Iodomctric Redox titration. 
Low field ac susceptibility measurements arc done with 



TABLE I: Structural parameters are determined from the Ri- 
etveld refinement of the powder XRD patterns for nanocrys- 
talline Pro.sSro.sMnOa. Here Ol is the apical and 02 is the 
equatorial oxygen in perovskite structure. Dxrd and Dtem 
are the average particle sizes calculated from the XRD and 
TEM measurements. Tb is the blocking temperature ob- 
tained from ac susceptibility. The extreme right column shows 
the percentage change (A%) in structural parameters for all 
the nanoparticles. 



Samples 


N600 


N650 


N700 


A% 


a (A) 


5.4270(7) 


5.4259(5) 


5.4255(5) 


-0.02 


c(A) 


7.7267(17) 


7.7226(12) 


7.7213(12) 


-0.06 


v (A 3 ) 


227.57(6) 


227.35(5) 


227.29(5) 


-0.12 


Mn-01 (A) 


1.9317 


1.9306 


1.9303 


-0.07 


Mn-Ol-Mn 


180° 


180° 


180° 


0.0 


Mn-02 (A) 


1.9337 


1.9278 


1.9303 


-0.25 


Mn-02-Mn 


165.7° 


168.6° 


167.1° 


+1.7 


Dxrd (nm) 


15.7 


17.3 


19.1 




Dtem (nm) 


15.7 


19.2 


26.6 




Tb (K) 


209.8 


240.3 


253.3 





home made ac susceptometerj* . DC magnetization have 
been measured with homemade vibrating sample magne- 
tometer (VSM)2I and commercial 14 Tcsla VSM (PPMS) 
made by Quantum Design. 



III. RESULTS AND DISCUSSIONS 
A. Superparamagnetic behavior of nanoparticles 

Fig. 2 presents the temperature variation of dc mag- 
netization (M) measured in 100 Oc following zero field 
cooled (ZFC) and field cooled warming (FCW) proto- 
cols. ZFC magnetization shows broad peak at character- 
istic temperature (Tb ) which is unlike its bulk compound 
where the sample exhibits clear PM-FM and then FM-AF 
transition on cooling from room temperature ! 14 ' 22 The 
observed peak which is typical feature of mctastable mag- 
netic systems like SPM particles or SG system, shifts to 
higher T with increasing particle size. It is clear in figure 
that ZFC and FCW branches of magnetization split at 
temperature (Ti rr ) higher than Tb, and Ti rr reduces with 
increasing the field (not shown). The higher value of Ti rr 
than Tb indicates the presence of larger size spin clusters 
which are ordered at higher temperature. Contrary to the 
bulk PSMO, the first order FM-AF transition at lower 
temperature as well as the associated hysteresis between 
FCW and the field cooled cooling (FCC) magnetization, 
down to 2K, is not evident in these nanoparticles. Similar 
phenomenon, however, is observed in various nanocrys- 
talline half doped manganites^— From the Curie- Weiss 
fitting of high-T magnetization data [(x = M/H)^ 1 vs 
T], the effective PM moment [n e ff) has been calculated 
in terms of Bohr magneton (/j,b) per formula unit (f.u.). 
The obtained values are 3.62, 3.70 and 4.38 ^s/f.u. for 
N600, N650 and N700 compound, respectively. It can 
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be noted that /Li e // for smaller particles are lower than 
the expected spin only value \p. e ff = g\/ S(S + 1), where 
S is the total spin] which is 4.38 /i B /f.u. for PSMO. 
It is noteworthy that with increasing particle size, ^eff 
approaches the expected value. 

To understand the low-T magnetic state in these 
nanoparticles we have collected high held (120 kOe) M 
vs H plots at 2 K (Fig. 3a). Unlike the bulk PSMO 
which at low T (within AF state) exhibits field induced 
jump in M(H) and attains the saturation moment (3.5 
/xb/i-u.) , 15 ' 23 the nature of M{H) plots for nanoparticles 
resemble qualitatively with FM systems. However, even 
in high field, magnetization does not saturate but show 
monotonous increase. The magnetic moments at 2 K and 
120 kOe (/it) are 1.38, 1.40 and 1.94 fj, B /lu. for N600, 
N650 and N700 sample respectively. These values are 
much lower than the saturation moment of bulk PSMO 
(3.5 //s/f.u.). Nonetheless, the estimated moments imply 
that the increased surface to volume ratio at smaller par- 
ticles creates higher amount of magnetically dead layer 
at the surface, thus it reduces the average moment sys- 
tematically. In Fig. 3b, M(H) data are plotted in the 
form of M 2 vs H/M which is known as Arrott ploti 2 ^ 
Positive intercept on the M 2 axis from the extrapolation 
of high field data in Arrott plot implies the presence of 
spontaneous magnetization in system. The spontaneous 
magnetization calculated from these intercepts (Fig. 3b) 
are 0.943, 1.003 and 1.367 £i B /f.u. for N600, N650 and 
N700 sample respectively. This clearly indicates that the 
ground state of these nanoparticles are FM contrary to 
their bulk counterpart even though there is no change in 
the structural symmetry. 

To confirm the magnetic state of these nanoparticles 
we have measured ac susceptibility (x) which probes the 
dynamics of the spin system in very low field. In general, 
magnetization can be expressed in terms of measuring 
field as: 

M = M + XiH + X2H 2 + xaH 3 + XiH A + .... (1) 

where Mo is the spontaneous magnetization, xi 1S the 
linear and X2, X3j Xi are the nonlinear susceptibilities. 
Fig. 4a presents temperature variation in real part of 
Xi (x?) measured in ac field 3 Oe and frequency 731 
Hz. The feature of Xi is similar to its dc magnetization 
counterpart presented in Fig. 2. The T B calculated from 
low field ac~x are mentioned in Table 1. 

It is rather nontrivial issue to distinguish between the 
SPM and SG as both share similar kind of experimental 
features. To do this, we have used nonlinear ac suscep- 
tibilities which are shown to be an useful experimental 
probes to separate out different magnetic states! 15 ' 25 " — 
In Fig. 4b, we have plotted temperature dependence of 
second order (X2) ac-x measured in ac field 3 Oe and fre- 
quency 731 Hz. Usually, X2 appears due to the presence 
of symmetry breaking field which originates either from 
the spontaneous magnetization or the superimposed dc 
magnetic fiekh 15 ' 26 ' 27 As the canonical SG does not hold 
spontaneous magnetization, so X2 should be absent with- 



out dc field. The presence of finite X2 in Fig. 4b confirms 
the presence of FM interaction within these nanoparticles 
and discards the possibility of canonical SG like phase. 
However, it is important to identify if the broad peak 
in ZFC magnetization and in x?" ca n be attributed to 
cluster-glass like phase or SPM. This is achieved through 
the measurement of 3rd order ac-x [X3(^)] which show 
broad negative peak around T B for all the samples (not 
shown). The SPM like feature in the present nanopar- 
ticles is confirmed from the noncritical behavior of the 
X3 around the peak temperature. It has been unambigu- 
ously shown that, for SG like system, X3 diverges as H 
— > around the transition temperature , 25 1 30 In Fig. 5, 
we have plotted the maximum value of X3 normalized by 
that at H = 20 Oe as a function of measuring field. As 
evident in figure, X3 is riot divergent as H — > 0. This 
analysis shows conclusively that PSMO nanoparticles do 
not show SG like coperative freezing but thermal block- 
ing of the magnetic entities. 

Further, the SPM behavior in these nanoparticles is 
also substantiated from the Wohlfarth's modeL 28 ' 31 The 
magnetization for an assembly of particles is given as: 

M = n(n)L((fi)H/k B T) (2) 

where n is the number of particles per unit volume, (//) 
is the average moment of magnetic particle, k B is the 
Boltzmann constant and L(x) is the Langevin function. 
After the expansion of L(x), xi and X3 above the block- 
ing temperature (T B ) can be expressed as: 

Xi =n(nf /'3k B T = P 1 /T (3) 

X3 = (n (M> /45)«m) /k B T f = P 3 /T 3 (4) 

Therefore, xi and X3 above T B for SPM will vary as T~ 4 
and T~ 3 , respectively. The main panels of Figs. 6a and 
6b show this variation of X\ and Xs following Eqs. 3 and 
4, respectively for N650 and N700 samples. In the inset 
of both the figures we show the same plot for N600 sam- 
ple. This further corroborates the earlier analysis that 
peak in xi arises due to the blocking of magnetic clus- 
ters in the present study. From the ratio of P3 and Pi, 
(/i) is estimated which are around 14.5 x 10 4 , 22.9 x 
10 4 and 26.1 x 10 4 fji B for N600, N650 and N700 sample, 
respectively. This large value of (/Li) is consistent with 
the SPM clusters as it consists of large number of spins, 
whereas for PM, (/Li) signifies only ionic moments and is 
limited to few fi B . Assuming that clusters are spherical, 
effective sizes of them are calculated from the values of 
(/z) derived from above analysis. These sizes arc found 
to be around 16.5, 19.2 and 20.9 ran for N600, N650 and 
N700, respectively. It is noteworthy that effective size of 
clusters calculated from the magnetic measurement are 
close to the same obtained from XRD and TEM mea- 
surements (Table 1). This consistency of particle size 
calculated from various methods proves the unambiguity 
in present analysis. 
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B. Interaction among the nanoparticles 

After having confirmed the SPM nature of these 
nanoparticles, we attempt to establish the existence of 
interparticlc interaction. As a first step, intcrparticle in- 
teraction is identified from the temperature dependence 
of coercive field (He) estimated from the M(H) plot for 
all the samples. For the noninteracting SPM particles, 
temperature variation in He is expressed asi22 



He — H, 



C,o 



(1 - T/T B ) 



1/2 



(5) 



where Hc,o is the value of He when T — > 0. Thus, the 
measured He as a function of T 1 / 2 would be linear for 
noninteracting particle systems. However, nonlinear be- 
havior of He vs T 1 / 2 in Fig. 7 arises due to the interac- 
tion among the particles. It is also evident in the figure 
that with increasing particle size nonlincarity of plot is 
reduced, thus probably indicate that the strength of in- 
teraction reduces with the increasing particle size. 

Now we attempt to reestablish the presence of inter- 
action among the particles from the frequency (/) de- 
pendent shift in Tb in ac-x(T) measurements. It can 
be mentioned that the variation of / in ac-x renders 
the variation in probe time rp (oc /~ 4 ) which allows to 
probe the relaxation of particles in different time win- 
dows. It is observed in Fig. 8 that for N600 sample 
decreases and Tb shifts toward higher T with increasing 
/. However, the shift of Tg with the / in present study 
is relatively small. It can be mentioned that simple, non- 
interacting SPM particles show large / dependence of 
Tb whereas for SG or interacting SPM particles, Tb is 
less / dependent^ The / dependence of Tb is normally 
quantified with the following empirical relation: 



AT fi 



T B Alog 10 (/) 



(6) 



where A is the difference in related parameters. The 
value of 4> for N600 sample is found to be around 0.002. 
The other two samples also show similar values of $. The 
experimentally determined values of $ for SPM particles 
are in the range of 0.1 - 0.13, whereas much lower values 
(0.005 - 0.05) are observed for canonical SG systems as 
well as interacting particle system . 34 ! 35 Therefore, based 
on the value of 3> it is often difficult to distinguish be- 
tween SPM and SG experimentally. However, the calcu- 
lated $ hints toward the presence of interaction among 
these particles. 

To unambiguously identify the presence of intcrparti- 
cle interaction, we follow a rigorous systematic method 
where the / dependence of Tb has been analyzed with 
different phenomenological models. According to Neel- 
Arrhennius law, relaxation time (r) for an assembly of 
noninteracting SPM particles behaves likej 3 ^ 



( E <* \ 
t exp I - — = I 



\k B TJ 



(7) 



E a is the anisotropy energy barrier equal to KV, where 
K is the anisotropy energy constant and V is the volume 
of the clusters. For the SPM relaxation, the prefactor 
t is in the range of 10 -8 - 10~ 13 sM Inset of Fig. 8 
shows the straight line fitting of Eq. 6 for N600 sample. 
Though fit is apparently good but the obtained values 
are quite unphysical (to ~ 10 -522 and E a /ks = 253615) 
which clearly indicates that dynamics of these nanopar- 
ticles can not be explained with noninteracting particle 
model. Thus, we have tried to analyze the / dependent 
shift of Tb using Vogel-Fulcher lawi2& 



t exp 



E a 



k B (T-T ) 



(8) 



The term T is the characteristic temperature (0 < T < 
Tb) which accounts the interaction among the particles. 
In this case, best fitting of data for N600 has been shown 
in Fig. 9a. The obtained fitted parameters are quite 
reasonable with the values E a /kB = 25.27 ± 3.34 K, To = 
205.15 K and t = 1.76 x 10" 8 s. Such higher value of t 
has been observed for interacting clusters^ This fitting 
confirms the interacting nature of these nanoparticles, 
and is in conformity with the recent observations like slow 
magnetic relaxation with logarithmic time dependence 
and memory effects^ 

It is rather significant and intriguing that the / de- 
pendent shift of Tb can be fitted to the scaling law as 
well, which is used to characterize the phase transition 
in SG, though the present system are proved to be al- 
most tailor-made small particles, which behave like SPM 
beyond doubt. The scaling hypothesis assumes that the 
relaxation time (t) is related to the correlation length 
(£) near to the transition temperature (T g ). As £ di- 
verges at T g , relaxation time obeys the following empiri- 
cal relation*^ 



t = T (T/T g - 1)- 



(9) 



where z is the dynamical scaling exponent and v is the 
critical exponent related to £. Fig. 9b shows the best 
fit of Eq. 8 for N600 sample and the obtained parame- 
ters are t = 3.39 x 10~ 14 s, T g = 206.44 K and zv = 
4.52. Usually, the exponent zv shows large variation for 
different kind of systems^ Although, the obtained zv is 
close to the theoretically predicted value (4) for 3D Ising 
model^S and the experimentally calculated value (5.5) 
for canonical SG CuMn 4 . 6 %i^ However, t shows large 
disagreement between Vogel-Fulcher law and scaling law 
fitting. Further, To obtained in later case is orders of 
magnitude less than that usually seen for SG (10 -11 - 
10 -12 s). This inconsistency of To with that for SG is 
rather significant, which discards the possibility of SG 
like freezing in these nanoparticles.— It is fact that the 
significant presence of intcrparticle interaction induces 
the collective behavior and glassy dynamics even in the 
tailor made SPM system. Since the response of such 
systems is similar to SG, they are termed as super spin 
glass (SSG) i 13 ' 41 " — Basically, SSG dynamics differ from 
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the SPM one in terms of critical slowing down of relax- 
ation time. Similar critically slow dynamics resembling 
SG like response is also shown by magnetic clusters form- 
ing in a matrix out of composition fluctuation or atomic 
segregation which is called cluster glass (CG) 4^ However, 
it remains an experimentally challenging task to discern 
between critically slow dynamics of SSG and CG with 
the slow but non-critical dynamics of SPM system. The 
high value of To and To obtained from the Vogcl-Fulcher 
law fitting unambiguously imply the presence of inter- 
particle interaction in present nanoparticles^ however, 
the obtained To from the dynamical scaling hypothesis 
fitting is a few orders lower than the value for SSG or 
CG (t = 1(T 9 - Kr 6 ).i^-i2 Nonetheless, the fitting 
of scaling law for our data is quite intriguing indicating 
some sort glassincss in the system, however, the overall 
results do not comply with the SSG behavior. This indi- 
cates that a careful exercise is needed to make a proper 
distinction between the SSG/CG like phases with that of 
SPM system from the dynamical behavior. 



C. Nature of interparticle interaction and its 
consequences on magnetic properties 

To identify the nature of the interparticle interaction 
and its consequences on the magnetic properties, we have 
measured the remnant magnetization following isother- 
mal remnant magnetization (IRM) and dc demagnetiza- 
tion (DCD) protocols^ For the IRM measurement, we 
cooled the sample in zero field from the room tempera- 
ture to 150 K and applied a field (H a ) isothermally for 
10 s. Then the H a is reduced to zero and remnant mag- 
netization is measured. We heated the sample back to 
room temperature and repeated this experiment for dif- 
ferent H a at the same temperature, each time increasing 
H a in steps of 50 Oe up to maximum of 1000 Oe. The 
measured IRM as function of H a for N600 sample is de- 
picted in Fig. 10a. In DCD measurement, the sample is 
zero field cooled from room temperature to 150 K, and 
magnetized to saturated state with applying +1000 Oe 
for 10 s. The applied positive field is made zero, and 
isothermally a negative field {—H a ) is applied for 10 s. 
Then, — H a is switched off and remnant magnetization is 
measured. We repeated this experiment similar to IRM 
measurement for different — H a up to -1000 Oe in step of 
-50 Oe. The measured DCD as function of \H a \ for N600 
sample is depicted in Fig. 10a. It can be mentioned that 
IRM and DCD originate from the virgin and saturated 
state, respectively. 

To characterize as well as quantify the magnetic in- 
teraction, the parameter SM is shown to be very useful 
which is defined as. 45 ' 46 

SM = m D cD - (1 - 2m mM ) (10) 

where m is the normalized remnant magnetization with 
respect to the remnant magnetization measured after the 



sample is magnetized to saturated state. For noninter- 
acting particles SM is zero4£ However, finite interactions 
among the particles lead to deviation of SM from zero, 
i.e., positive SM is due to the interactions which favor 
magnetization whereas negative SM arises from the in- 
teractions which cause demagnetization in the system4£ 
Usually, the negative SM implies the dipolar interaction 
and the field where minimum in SM occurs quantifies 
the strength of interaction^ We have plotted SM as a 
function of \H a \ for N600 sample in Fig. 10b. As evident 
in figure, SM is negative in low H a , showing minimum 
around 150 Oe. In higher H a around 522 Oe, SM shows 
crossover to positive value. We can infer from this results 
that interparticle interaction is of dipolar type, having 
magnitude roughly of 150 Oe. At high field, interparti- 
cle interaction is dominated by the applied field yielding 
positive SM. For the same measurements on N700 com- 
pound, SM shows negative minimum around 100 Oe, and 
crossover to positive value around 220 Oe. These results 
indicate that with the increasing particle size, strength of 
dipolar interaction reduces. The dipolar interaction be- 
tween two magnetic dipoles pf, and /Tj can be estimated 
as)2£ 



Ed = yti-fo - 3 [&-hj) (fij-hj)] (11) 

ij 

where is the center to center distance between the 
dipoles. Eq. 11 gives the form of long-range dipolar 
interaction energy between the two magnetic diploes of 
moment (/x) separated by distance . For a collection of 
small particles, it has to be summed over all the particles. 
The present system may be considered as an assembly of 
spherical nanoparticles of radius R, placed adjacent to 
each other in a regular array. In this case, the sepa- 
ration between the adjacent magnetic moments, rij, is 
about 2R. If we consider the magnetic moment (fi) of 
each particle is roughly proportional to its volume, then 
according to Eq. 11, dipolar interaction energy (Ed) be- 
tween the two adjacent particles varies approximately as 
R 3 . Since for a collection of particles within a given vol- 
ume, the number of particles varies as inverse of R 3 , the 
total interaction energy which is the summation over the 
entire individual E^s, becomes almost independent of the 
particle size. However, in real situation, as the average 
particle size decreases, the packing fraction increases in 
such densely packed assembly of nanoparticles having a 
distribution in particle size. Thus decrease in the size 
will result in the increase in the number of particles in 
the same volume, much more than what is expected for 
a regular array where it increases as 1/R 3 . Moreover, for 
such dense packing, reduction in the size of the particles 
will also reduce the average distance between them, r^, 
faster than a linear function of R. Both this factors will 
contribute to increase the total dipolar interaction energy 
in the system with the decrease in particle size. The evi- 
dence of such increase in dipolar interaction energy with 
the decrease in the particle size in the present system can 
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be found in the Fig. 7 and the discussions of sections B 
and C. 

The dipolar interaction has significant influences on 
the magnetic properties as it modifies the otherwise uni- 
axial anisotropy energy barrier of individual nanoparti- 
cle. The most notable effects of this interaction are the 
evolution of glassy dynamics ji 3 - modification in Tb ) 13 i 48 
etc. Here we have looked into the effect of dipolar in- 
teraction on the field dependence of Tb by measuring 
ZFC magnetization in different fields. In the experimen- 
tal scenario, Tb corresponds to the temperature where rp 
equals to the average r of the system. Since with the vari- 
ation of applied field Tb varies, thus one can construct a 
line of constant r on the plane spanned by H and T. The 
existence of such line is predicted in the mean-field theo- 
retical model for SG. To be specific, this model assumes 
phase transition in SG, and considers that large field will 
destroy the frozen spin stated Therefore, in presence of 
field, the critical lines are predicted on the H-T plane 
which mark the phase transition. The first one is the 
de Almeida-Thouless (AT) line occurring in anisotropic 
Ising SG and behaves as Tb(H) oc H 2 / 3 . The second 
one is the Gabay- Toulouse (GT) line valid for isotropic 
Heisenberg case, and shows functional form Tb(H) oc 
j{2 36 However, later numerical calculation showed that 
such lines are not unique for SG, and can even exist in 
case of relaxation of (interactive) SPM particles with the 
crossover from low-H GT to high-iJ AT behavior^ Our 
results arc plotted in Fig. 11 where we show that vari- 
ation of Tb with field follows AT line in low field for 
all the samples. However, with the increase in field we 
find the deviation from AT line as marked by vertical 
arrows in the figure. Remarkably, in high field we find 
the variation of Tb (H) agrees with the GT line. We have 
observed the similar behavior for all the samples and in 
inset of Fig. 11 we have plotted such high field behavior 
for N600 sample. The fields (H cr ) where we find crossover 
from AT to GT behavior are approximately 250, 200 and 
100 Oe for N600, N650 and N700 sample, respectively. 
It is rather significant that the crossover fields from AT 
to GT like behavior are close to the dipolar interaction 
fields found from the remnant magnetization measure- 
ments mentioned above (Fig. 10). This crossover in crit- 
ical lines around the dipolar interaction field for these 
nanoparticles is quite noteworthy. 

The above experimental results in Fig. 11 has the 
similarity with the theoretical calculation by Kotliar and 
Sompolinsky 5 ^ who predicted that in presence of unidi- 
rectional random anisotropy (i.e., Dzyaloshinsky-Moriya 
type interaction) the critical behavior for SG in fields 
lower than the anisotropy is close to Ising-likc follow- 
ing AT line, and crosses over to Heisenberg behavior in 
high fields. In fact, similar anisotropy induced crossover 
in critical lines is experimentally observed in Au doped 
classical CuMn SG system as well as in superconduct- 
ing vortex-glass system.] 51 ! 52 These results allow us to 
draw an analogy with the present PSMO nanoparticles 
where the interparticlc interaction acts as an unidirec- 



tional anisotropy leading to crossover from AT to GT 
like behavior at field around the interaction field. Indeed, 
dipolar interaction (Eq. 11) is anisotropic which intro- 
duces angular dependence of spin ordering in system^ 



To confirm the unidirectional nature of the anisotropy 
arising out of interparticle interaction we have measured 
field cooled (FC) M vs H loop for all the nanoparti- 
cles. In presence of unidirectional anisotropy, FC hys- 
teresis loop is shifted along the field axis, generally, in 
opposite direction to the cooling field. This loop shift 
is commonly known as 'Exchange Bias (EB)' which is 
characterized by field, He^ The samples under study 
have been cooled in +10 kOe from room temperature to 
10 K, and after proper thermal stabilization hysteresis 
loops have been recorded (Fig. 12). We have ensured 
that the artifacts due to minor hysteresis loop arc elimi- 
nated. It is evident in figure that collected M vs H loop 
is shifted toward the negative field axis. It is worth men- 
tioning that we find similar opposite shift when samples 
are cooled in negative field. We calculate He [= - [h + 
- h-)/2, where h_ and h+ are the point of intersection 
on the field axis at decreasing (— ) and increasing (+) 
fields cycles] as 155.3, 132.4 and 105.7 Oe at 10 K for 
N600, N650 and N700 sample, respectively Convention- 
ally, EB is believed to be directly associated with the 
interface of FM and AF components. However, FM in 
contact with SG or ferrimagnets are also observed to give 
rise EB properties. 53 ' 54 For the present PSMO nanopar- 
ticles, we believe that the observed EB is guided by the 
interparticle interaction rather than the FM/AF inter- 
face mechanism. Our conclusion is based on the follow- 
ing points: (i) Fig. 3 conclusively shows the FM nature 
of low-T magnetic state for all the nanoparticles. (ii) 
Even if, EB arises due to the presence of residual AF 
components then He would vanish above the FM-AF 
phase transition temperature Tjy. We find finite He at 
or above 150 K whereas Tjy during cooling in same field 
occurs around 84 K for bulk PSMOi22 (iii) The resid- 
ual AF components arc supposed to increase with the 
increasing particle size as the decrease in particle size 
does not favor AF ordering therefore it will result 
in more FM/AF interfaces. But, our estimated He de- 
creases with increasing particle size, (iv) The variation 
of He follows the similar trend of interparticle interac- 
tion, i.e., decreases with increasing particle size. These 
results straightforwardly imply that EB in these PSMO 
nanoparticles originate due to the mechanism guided by 
the interparticle interaction. Furthermore, these confirm 
the unidirectional-anisotropic nature of this interparticle 
interaction which leads to crossover in critical lines on 
the H-T plane. Details of EB in these nanoparticles will 
be published elsewhere. Nonetheless, this interparticle 
interaction induced crossover in critical lines with fields 
in these compounds is quite intriguing and requires fur- 
ther studies involving both experimental and theoretical 
endeavor to comprehend. 
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IV. CONCLUSION 

In conclusion, we show that nanoparticles of 
Pr .5Sr .5MnO3 of three different average sizes have FM 
ground state contrary to their bulk counterpart in spite 
of having same crystallographic symmetry. The low tem- 
perature FM-AF transition which is a marked feature 
of bulk compound is significantly absent in nanoparti- 
cles. Detailed linear as well as nonlinear ac suscepti- 
bilities coupled with dc magnetization confirm the SPM 
nature of these nanoparticles. Presence of interactions 
among these particles has been corroborated from the 
analysis of frequency dependent peak in ac susceptibil- 
ity and temperature dependence of coercive field. The 
nature of this interaction is identified to be of dipolar 
type and the strength of the interaction is found to de- 
creases with the increase in particle size. This interpar- 
ticle interaction gives rise to some kind of glassiness in 
the magnetic response even in these SPM system. The ef- 
fect of this dipolar interaction on the magnetic properties 



is clearly evident as the systems exhibit crossover from 
AT to GT like critical lines with increasing field above 
their respective interaction field, which is rather intrigu- 
ing. Following the theoretical prediction, we believe that 
this crossover phenomenon is induced by the presence of 
unidirectional anisotropy which arises from the interpar- 
ticle interaction and the presence of exchange bias effect 
in these samples confirms it. 
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FIG. 1: (Color online) (a) Rietveld refinement profile of powder XRD pattern of N650 sample at room temperature. The 
open circles and solid continuous lines represent the observed and calculated pattern respectively. The difference plot is shown 
at the bottom of the figure, (b) The bright field TEM image of N650 sample collected at room temperature, (c) Histogram 
obtained from several TEM images shows particle size distribution for N650 sample, (d) Room temperature selected area 
electron diffraction (SAED) pattern of N650 sample showing good crystalline nature of nanoparticles. 
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FIG. 3: (Color online) (a) Magnetization vs field data are shown for all the nanoparticles at 2 K. (b) Arrott plot (M 2 vs H/M) 
of the isotherms in (a) are shown. Lines are due to straight-line fitting of the plots in high fields. The data for N650 sample 
are vertically shifted by 200 for clarity. 




FIG. 4: (Color online) (a) Real part of 1st order ac susceptibility measured in field 3 Oe and frequency 731 Hz are plotted 
as a function of temperature for all the nanoparticles. (b) Real part of 2nd order susceptibility is plotted as a function of 
temperature for all the nanoparticles in same field and frequency. 
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FIG. 5: (Color online) The maximum value of \3 normalized by that at 20 Oe has been plotted as a function of applied ac 
field for all the nanoparticles which clearly shows the non-diverging behavior of X3 as H — > 0. 
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FIG. 6: (Color online) Temperature variation in (a) 1st order and (b) 2nd order ac susceptibility is plotted above the blocking 
temperature for N650 and N700 samples. Straight lines are the (a) T _1 fit to \i (Eq. 3) and (b) T~ 3 fit to \3 (Eq. 4) data. 
Inset shows the same plotting for N600 sample in the respective graphs. 
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FIG. 8: (Color online) Temperature dependence of real part of xi at various frequencies is plotted for N600 compound. Arrow 
indicates peak shift with the increasing frequencies. Inset shows Neel-Arrhennius law (Eq. 6) fitting of frequency dependent 
peak temperature (defined in text). 
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FIG. 9: (Color online) The best fit of the relaxation times (r) to (a) the Vogel-Fulcher law (Eq. 7) and (b) the Scaling law 
(Eq. 8) for N600 sample. 
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FIG. 10: (a) Remnant magnetization measured following IRM and DCD protocols (denned in text) are plotted as a function 
of applied field for N600 sample. Bottom panel shows calculated SM values following Eq. 10 for the same sample. 
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FIG. 11: Color online) The temperatures where the peak occurs in ZFC magnetization are plotted as a function of applied field 
for all the samples. Straight line is fit to AT line (defined in text). The vertical arrows mark the field above which systems 
deviate from AT line behavior. Inset: Field dependence of peak temperature has been fitted to GT line for N600 sample 
(defined in text) 
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FIG. 12: (Color online) Field cooled M vs H loops collected after cooling the samples in 10 kOe from room temperature to 10 
K. The recorded loops exhibit shift along the field axis toward negative direction. 



